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(Dated: January 5, 2007) 

We have performed the first direct measurement of the time integrated flavor untagged charge 
asymmetry in semileptonic B® decays, Ag2 nt , by comparing the decay rate of B^ — ► ji + DjvX, 
where D^ — > (f>-K~ and <j> — > K + K~ , with the charge-conjugate B° a decay rate. This sample was 
selected from 1.3 fb _1 of data collected by the DO experiment in Run II of the Fermilab Tevatron 
collider. We obtain A"gl nt = [1.23 ± 0.97 (stat) ± 0.17 (syst)] x 10~ 2 . Assuming that Am s /f s > 1, 



4 



this result can be translated into a measurement on the CP-violating phase in B® mixing: Ar s /Am s 
tan 0s = [2.45 ± 1.93 (stat) ± 0.35 (syst)] x 10" 2 . 



PACS numbers: 12.15.Hh, 13.20.He, 14.40.Nd 



This letter presents the first measurement of a time 



integrated flavor untagged charge asymmetry Ag L in 
semileptonic B® decays. This asymmetry is defined as: 



A 



N(vl+D-)-N(»-D+) 
N(n+D7)+N(n-D+)' 
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where N(fj, ± Df) is the number of decays B° s — > 
H^DJvX integrated over the B° s lifetime. This asym- 
metry is called untagged because the initial flavor of the 
Bg meson is not determined. A s ^ nt is related to CP vio- 
lation in Bg mixing [l[ and can be expressed through the 
parameters of the B® mass matrix as 0]: 
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where Ar s (Am s ) is the width (mass) difference between 
the mass eigenstates in the B® system, x s = Am s /F s , 
y s = Ar s /(2f s ) where f s is the average width in the B® 
system, and <p s is a CP-violating phase. The standard 
model (SM) predicts a very small value for this asymme- 
try 2 x A S s U L nt = a s SL = (0.21 ± 0.06) x 10" 4 while the 
contribution of new physics can significantly modify this 
prediction 0, HJ . 

This measurement was performed using a large sample 
of semileptonic B® decays collected by the DO experi- 
ment at the Fermilab Tevatron collider in pp collisions at 
y/s = 1.96 TeV and follows closely the procedure used 
in the estimate of the dimuon asymmetry described in 
Ref. Q. The data correspond to an integrated luminos- 
ity of approximately 1.3 fb _1 . The DO detector is de- 
scribed in detail elsewhere 0]. The detector components 
most important to this analysis are the central tracking 
and muon systems. The central tracking system con- 
sists of a silicon microstrip tracker (SMT) and a central 
fiber tracker (CFT), both located within a 2 T supercon- 
ducting solenoidal magnet, with designs optimized for 
tracking and vertexing for pseudorapidities of \r)\ < 3 
and \r]\ < 2.5, respectively. The outer muon system, 
with coverage for \r/\ < 2, consists of a layer of track- 
ing detectors and scintillation trigger counters in front of 
1.8 T iron toroids, followed by two similar layers after the 
toroids @ . The polarities of the solenoid and toroids are 
reversed regularly during data taking, so that the four 
solenoid-toroid polarity combinations are exposed to ap- 
proximately the same integrated luminosity. The direct 
and reverse magnetic fields in the magnet were measured 
to be equal to within 0.1%. The reversal of magnet po- 
larities is essential to reduce the detector-related system- 
atics in asymmetry measurements and is fully exploited 
in this analysis. 



The asymmetry A S g™ 1 was measured using the decay 
B° -> nD s vX with D s -> <pir, <p -> K + K~. The se- 
lection of this final state is described in detail in Ref. 
Q. No explicit trigger requirement was made, although 
most of the sample was collected with single-muon trig- 
gers. Muons were required to have transverse momen- 
tum pt(p) > 2 GeV/c and momentum p(/j,) > 3 GeV/c, 
to have hits in both the CFT and SMT, and to have 
measurements in at least two layers of the muon sys- 
tem. All reconstructed charged particles in the event 
were clustered into jets [l(J, and the D s candidate was 
reconstructed from three tracks found in the same jet 
as the reconstructed muon. Oppositely charged parti- 
cles with pt > 0.7 GeV/c were assigned the kaon mass 
and were required to have an invariant mass 1.004 < 
M(K+K-) < 1.034 GeV/c 2 , consistent with that of a 
meson. The third track was required to have px > 0.5 
GeV/c, a charge opposite to that of the muon charge, 
and was assigned the pion mass. The three tracks were 
required to have hits in the CFT and SMT and to form 
a common D s ve rtex using the algorithm described in 
detail in Ref. ll| . To reduce combinatorial background, 
the D s vertex was required to have a positive displace- 
ment in the transverse plane, relative to the pp collision 
point (or primary vertex), with at least 4<j significance. 
The cosine of the angle between the D s momentum and 
the direction from the primary vertex to the D s vertex 
was required to be greater than 0.9. The trajectories of 
the muon and D s candidates were required to originate 
from a common B® vertex, and the (/xZ) s ) system was 
required to have an invariant mass between 2.6 and 5.4 
GeV/c 2 . 

To further improve the B° signal selection, a likeli- 
hood ratio method [13] was utilized. Using background 
sidebands (B) and sideband-subtracted signal (S) distri- 
butions in the data, probability distributions were found 
for a number of discriminating variables. These variables 
were the angle between the D s and K momenta in the 
K + K~ center-of-mass frame, the isolation of the (fiD s ) 
system, the \ 2 of the D s vertex, the invariant masses 
M{nD s ) and M(K + I<-), and p T {K+ K~). The isolation 
was defined as the ratio of the sum of the momentum of 
the tracks used to reconstruct the signal divided by the 
total momentum of the tracks contained within a cone 
with \J At; 2 + A<p 2 < 0.5 centered on the direction of the 
fj,D B system. The final requirement on the combined se- 
lection likelihood ratio variable was chosen to maximize 
the predicted ratio S/y/S + B. 

For this analysis we required the B® vertex to have a 
positive displacement from the primary vertex to sup- 
press the combinatoric background from the process 
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The numbers of events riq 7 (Z) s 
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FIG. 1: The invariant mass distribution M(07r) for the se- 
lected Bs candidates. The curve shows the result of fit by a 
function described in the text. 



cc(bb) — ► [iD s vX with the D s originating from a b or 
c quark, and the muon arising from another quark. The 
invariant mass distribution M (4>ir) for the selected events 
is shown in Fig. [TJ The low and high peaks correspond re- 
spectively to (/xD), mostly due to B°, and (/j,D s ), mostly 
due to B®. The curve represents a fit to the Al(4>ir) 
spectrum. A single Gaussian was sufficient to describe 
the D — > 4>ir decay, a double Gaussian to describe the 
D s — > (f>TT decay, and the background was modeled by an 
exponential. The total number of events passing all cuts 
in the D s mass peak is 27,300 ± 300 (stat). 



To measure A 



SL 



both physics and detector effects 
contributing to the possible imbalance of events with pos- 
itively and negatively charged muons must be taken into 
account. One physics source of asymmetry is CP vio- 
lation in scmilcptonic B decays. In addition, forward- 
backward charge asymmetry of events produced in the 
proton- antiproton collisions can also be present. Detec- 
tor effects can give rise to an artificial asymmetry if, e.g., 
the reconstruction efficiencies of positively and negatively 
charged particles are different. However, a positively 
charged particle produces the same track as a negatively 
charged particle in the detector with reversed magnet 
polarity. Therefore, almost all detector effects can be 
canceled provided the fractions of events with opposite 
magnet polarities are approximately the same. This is 
the case in this analysis, where the exposures are the 
same within 1%. 

According to the method described in Ref. 0, the 
event sample was divided into eight subsamples corre- 
sponding to all possible combinations of the toroid polar- 
ity = ±1, the sign of the pseudorapidity of the (n(j>n) 
system [13[ 7 = ±1, and the sign of the muon charge 
q = ±1. The number of (/J,D S ) events in each subsample 
was obtained by a fit to the mass distribution M((j)n) us- 
ing the same function as for the whole sample. For the 
cross-check we also extracted the numbers of (/J.D) and 



(£>)] in the 

D s [D] mass peak and in the background n^ 7 (bkg) for eight 
subsamples. 



Subsample: 






nf (bkg) 


13-yq 


(events) 


(events) 


(events) 


+ + + 


3,216 ±76 


907 ±55 


9,797 ±124 


H h 


3,586 ±79 


965 ± 56 


10,387 ±127 


+ + - 


3,391 ±78 


1,037 ±57 


10,390 ±127 


H 


3,225 ± 76 


963 ± 55 


9,832 ± 124 


- + + 


3,616 ±80 


1,003 ±57 


10,508 ±128 


h 


3,370 ± 77 


801 ± 54 


9,987 ±125 


- + - 


3,353 ±77 


831 ±55 


10,215 ± 125 




3,532 ±79 


1,116±59 


10,701 ±129 



background events from the fit. The widths and positions 
of the {nD s ) and (fJ,D) peaks, the relative fractions of the 
two Gaussians describing the (iiD s ) peak, as well as the 
background slope were fixed to the values obtained from 
the fit to the total M(07r) distribution. The numbers 
of (/J,D S ) and (fJ.D) events, n^{D s ) and n^ 7 (D), along 
with the number of the background events in the fitting 
range 1.75 — 2.30 GeV/c 2 , r^ 7 (bkg), for each subsample 
is given in Table fl] 

The fitted numbers of (fJ.D s ) [(pD), background] 
events were used to disentangle the physics asymme- 
tries and the detector effects. The n^ 7 can be expressed 
through the physics and the detector asymmetries as fol- 
lows I6|: 



4 v 



qA)(l + qjA {b ){l + jA Aet ) 
x (1 + q (3 1 A xo )(l + q(3A qP )(1 + (3jA 01 ) . 



(3) 



Here N is the total number of (nD s ) [(fiD), background] 
events; e 13 is the fraction of integrated luminosity with 
toroid polarity /3 (e + ±e~ = 1); A is the integrated charge 
asymmetry to be measured; Ar, is the forward-backward 
asymmetry; A^ct is the detector asymmetry for particles 
emitted in the forward and backward direction; A m is 
the range-out asymmetry that accounts for the change 
in acceptance of muons which bend towards the beam 
line and those which bend away from the beam line; A q p 
is the detector asymmetry which accounts for the change 
in the muon reconstruction efficiency when the toroid 
polarity is reversed; A^ accounts for any detector re- 
lated forward-backward asymmetries that remain after 
the toroid polarity flip. 

Since the system (|3]) contains eight equations, all six 
asymmetries together with N and e + can be extracted 
for each of the three types of the events. Results are pre- 
sented in Table [TTI separately for (nD s ) and (/i-D) events 
and the background. The physics asymmetries A and Afb 
for background events are consistent with zero. This is 
an important test for this method, since the precision of 
the asymmetry measurement for the background events 
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TABLE II: The physics and detector asymmetries for (fj,D s ), 
(l-iD) and background events. Uncertainties are statistical. 







(/iD) 


Background 


N 


27,289 ± 220 


7,623 ±162 


81,817 ±357 


6+ 


0.492 ±0.004 


0.510±0.011 


0.494 ±0.002 


A 


0.0102 ±0.0081 


-0.0345 ±0.0211 


-0.0056 ±0.0045 


An, 


-0.0046 ±0.0081 


0.0480 ±0.0210 


-0.0020 ±0.0043 


^4det 


-0.0051 ±0.0081 


-0.0072 ±0.0212 


0.0001 ±0.0044 


A TO 


-0.0352 ±0.0081 


-0.0819 ±0.0209 


-0.0263 ±0.0044 


A 0J 


-0.0097 ±0.0081 


0.0104 ±0.0213 


-0.0010 ±0.0044 


Kb 


0.0030 ±0.0081 


0.0014 ±0.0212 


0.0046 ± 0.0044 



is much higher than that of the signal due to the larger 
statistics. The largest detector asymmetry for all three 
types of the events is the range-out asymmetry. 

It can be seen from (j3|) that if e + = er = 1/2, and the 
asymmetries A, Afb, Adet> A IO , A q p, Ap^ are small, each 
of them can be obtained independently by the appropri- 
ate division of the entire sample of events into two parts. 
For example, the asymmetry A can be obtained be di- 
viding the sample according to the charge of muon. For 
such a division, and neglecting the second order terms, 
we obtain: 



.4 



(4) 



,07 



qA). 



This observation explains in particular the similar values 
of statistical uncertainties for all asymmetries in Table 

IE 

The resulting charge asymmetry of (fJ,D s ) events is 
A = 0.0102 ± 0.0081 (stat). It is related to A S g U L nt via 
A = f s ■ A s g1 nt + f d ■ A s '2 nt , where f s (f d ) is the fraction 
of B^(B^) — * fiD s i/X decays in the (//-D s ) sample. A s '2 nt 
may arise only from B® — > DD S decay, the fraction of 
which in the (nD s ) sample was found to be small, at the 
level of (4 ± 1)%. Additionally, the value of A s '2 nt is 
strongly constrained experimentally 3, HI to be close 
to zero. Therefore the time-integrated Ag^ nt compo- 
nent can be neglected. The fraction of B° s decays, f s , 
was determined as follows. The decays B® — > \xD s vX 
and B° s — > tD s vX — > [iD s vX were considered as a 
signal. The decays B® — > D S D S X with D s — > [ivX 
are not flavor-specific and hence were considered as a 
background. The decays B a d — ► DD S X were also in- 
cluded in the background. In addition, the process 
cc(bb) — ► fiD s vX was taken into account. This back- 
ground produces a pseudovertex which peaks around the 
primary interaction point. It is reduced by approximately 
50% by requiring a positive displacement of the {^D s ) 
vertex. 



All processes were simulated using the evtgen [16[ 
generator interfaced to PYTHIA 17[ and followed by full 
modeling of the detector response using GEANT [18| and 
event reconstruction. The branching fractions of B d de- 
cays were taken from Ref. [l[, while the contribution 
of the process cc(bb) — + [iD s vX was measured directly 
in our data to be (5.9 ± 1.7)%. With these assump- 
tions, (83.2 ± 3.3)% of the selected sample of {^D s ) 
events is composed of semileptonic B° s decays. The un- 
certainty on this value comes from the uncertainties on 
the branching ratios of the contributing B decays and 
the uncertainty on the fraction of the cc(bb) — > fj,D s vX 
process in the sample. Taking into account the sample 
composition, the measured integrated charge asymme- 
try of semileptonic B® decay is found to be Ag^ nt = 
[1.23 ±0.97 (stat)] x 10~ 2 . 

The following sources of systematic uncertainty were 
considered. The final state includes a K + K~ pair. 
Therefore, the charge asymmetry of K meson reconstruc- 
tion, which arises due to the different interaction cross 
sections of K + and K~ in the detector material, does not 
contribute to the measured A 



SL 



The charge asymme- 
try of pion reconstruction, however, can contribute. The 
nd interaction cross sections for positive and negative pi- 
ons differ by (1.3 ± 0.3)% in the range 1 - 2 GeV/c 
Taking into account the amount of material which a pion 
crosses in the detector, the induced asymmetry due to 
pion reconstruction was estimated to be 2 x 10~ 4 . This 
value was included in the systematic uncertainty. 

The uncertainty in the fraction of B° s signal in the 
(/i£> s ) sample produces a systematic uncertainty of 1 x 
10~ 3 . This uncertainty also includes a possible residual 
variation of the signal fraction between subsamples. 

The uncertainty due to the fitting procedure was es- 
timated by varying the masses and widths of the peaks, 
and the slope of the background by one standard devi- 
ation. The fitting procedure was also repeated with a 
single Gaussian describing the D s peak and with a dif- 
ferent fitting range. The resulting change of Ag^ nt did 
not exceed 0.14 x 10~ 2 which was used as an estimate of 
the systematic uncertainty from this source. 

The reconstruction efficiency varies with the de- 
cay length due to the applied requirements. We verified 
that this variation does not bias the result for j4c'^™ and 
the relation @. In addition, any possible contribution 
of the B® charge asymmetry to the measured value was 
estimated to be negligible. 

Adding all contributions into the systematic uncer- 
tainty in quadratures, we obtain the resulting value of 
the time-integrated untagged charge asymmetry: 



.4 



s,unt 
SL 



[1.23 ±0.97 (stat) ± 0.17 (syst)] x 10" 2 .(5) 



This is the first direct measurement of A S gj™ . It can be 
seen that the statistical uncertainty dominates and will 
be improved in the future with the increase of statistics 
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and addition of new decay modes. Using Eq. ((2|) and 
assuming that Am s /r s 3> 1, we obtain: 



tan( 



[2.45 ± 1.93 (stat) ± 0.35 (syst)] x 10" 



(0) 



This result, together with the measurements of AT S 
21 1 and Am s [ij, [22| , provides a constraint on the CP- 
violating phase 4> s . 
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